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THE CALCULATION OF CERTAIN gTATIC AEROELASTIC 
PHENOMENA OF WINGS WITH TIP TANKS OR 
BOCM-MOUNTED LIFTING SURFACES 
By Franklin W. Dlederlch and Kertneth A. Foss 

StMMARY 


The matrix- integration method of MCA Rep. 1000 for calculating 
static aeroelastic phenomena is extended to the case of a wing with 
concentrated aerodynamic forces at the tip due to tip tanks or hoom- 
mo\mted lifting surfaces. A simplified method of calculation which is 
based on the concept of the semirigid wing and which utilizes the pre- 
sumably known aeroelastic characteristics of the wing alone is presented 
for cases in which the aerodynamic interaction between the concentrated 
force and the remainder of the wing can be neglected. The modified 
matrix-integration method has been used to calculate some static aero- 
elastic characteristics of an unswept wing with a tip tank, and both 
methods have been used to calciilate the characteristics of a 45° swept- 
back wing with several boom-mounted lifting- surface configurations. 

The results of these calc\ilations show that the presence of a tip 
tank on an unswept wing tends to deteriorate its static aeroelastic 
characteristics and that a lifting surface geared to the aileron and 
mounted on a boom ahead of the tip of a sweptback wing may improve the 
static aeroelastic characteristics of the wing to a sufficient extent 
to warrant consideration of such a vane as a device for relieving 
adverse aeroelastic effects. 


INTRODUCTION 


The matrix- Integration method of references 1 and 2 for calculating 
static aeroelastic effects of swept wings of arbitrary stiffness implies, 
as do most other methods of static aeroelastic analysis which treat the 
wing essentially as a simple beam, that the vertical shear, moment, and 
torque at the tip are zero. The presence of concentrated forces and 
moments at the wing tip violates this assumption to the extent that 
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they are discontinuous at the tip. One case in which such forces are 
of interest is a wing with a tip tank. Another case is a wing with a 
hoom-mounted lifting surface. 

This combination may he of interest because a surface mounted on a 
boom ahead of the tip of a sweptback wing Introduces large twisting 
moments and, if the surface is mounted such that its angle of attack is 
the same as that of the wing tip, it causes twisting of the structure 
in a direction such as to oppose the effect of the bending deformations; 
hence, by reducing the net change of angle of attack due to wing defor- 
mation the vane tends to reduce the shift of the aerodynamic center due 
to aeroelastic action. Since aerodynamic forces due to an aileron 
deflection cause twisting and bending deformations both of which give 
rise to aerodynamic forces which tend to oppose those due to the aileron 
deflection. Increasing the twisting deformation by means of a boom- 
mounted surface only tends to aggravate the loss of lateral control due 
to aileron deflection. However, if the siirface is geared to the aileron, 
so that it pitches up when the aileron is deflected downward, it tends 
to reduce the amoiint of lateral control lost because of aeroelastic 
action. Furthermore, it may increase the lateral-control power sub- 
stantially under certain conditions when there is no aeroelastic action, 
as, for instance, when the aileron is relatively ineffective because 
of boundary- layer accumulation or because of shock on the wing ahead of 
the aileron. Consequently, a boom-mounted geared lifting surface 
appears to warrant consideration as a device for alleviating adverse 
ae 2 X>elastic effects. 

For these reasons the method of reference 1 is extended to the 
case of concentrated forces at the wing tip in the present paper. In 
this modified method, most of the matrices used in the analysis of the 
wing alone by the method of reference 1 can also be used in the calcu- 
lations for the wing with the concentrated force at the tip. If 
aerodynamic-induction effects between the wing proper and the body 
producing the concentrated aerodynamic force under consideration are 
neglected, a simpler method may be used to calculate the desired aero- 
elastic effects. Such a method is also described in this paper; it 
consists of correcting the presumably known aeroelastic effects of the 
wing alone for the presence of the concentrated force in a manner sug- 
gested by the semirigid-wing concept. 

In order to illustrate the results obtainable by these methods, 
calciilations have been made for an unswept wing with and without a tip 
tank and for a 45° sweptback wing with and without several boom-mounted 
lifting- siirf ace configurations. In the case of the sweptback wing, 
calculations have been made both by the matrix-integration and the 
simplified methods with substantially identical results. The results 
of the calculations are discussed and certain conclusions are drawn; a 
knowledge of the method of analysis is not required for an understanding 
of this discussion. 
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SYMBOLS 


A 


aspect ratio 


a location of local aerodynamic center rearward of leading 

edge, fraction of chord 

a location of wing aerodynamic center rearward of leading edge 

of mean aerodynamic chord, fraction of mean aerodynamic 
chord 


h wing span 

h ' wing span less fuselage width 

C _ , constants defined by equations ( 58 ) to ( 61 ) 

^'2 3 h constants defined by equations ( 69 ) to (71) 

Cj^ rigid-wing lift-curve slope per radian 

cc 

> 

Cj coefficient of dancing in roll (rolling-moment coefficient 

for linear antisymmetric twist of 1 radian at wing tip) 

Cu]^ pitching-moment-curve slope per radian 

c chord parallel to free stream 

"c average cho 2 rd (s/b) 

Cl section lift coefficient (z/qc) 

d distance parallel to free stream between center of pressure 

of boom-mo\mted lifting siorface and elastic axis 

El bending stiffness 

e local position of elastic axis rearward of leading edge, 

fraction of chord 

ei distance along chord from elastic axis to section aerodynamic 

center, fraction of chord (see fig. l) 

e 2 distance along chord from elastic axis to center of pressure 

due to aileron deflection, fraction of chord ( see fig. l) 
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GJ 

8>S' 

K 

% 

%,R,v 

k 


torsion stiffness 

factors defined "by equations (22) and (33) 
gear ratio between boom motion and aileron motion 
spring constant of boom 

coefficients defined by equations ( kO) , ( 4l ) , and ( 66) 

b V2 tan /A 

cos^Ay 

coefficients defined by eqixatlons (68) to (71) 


dimensionless parameter 



L lift 

1 section lift 

M accumulated bending moment (about an axis parallel to free 

stream, unless specified otherwise) 

Mq free -stream Mach number 

P concentrated normal force 

q dynamic pressure 

q* dimensionless dynamic pressure 

"q dimensionless dynamic pressure 

T accumvilated torsion moment (about an axis perpendicular to 

plane of symmetry, unless specified otherwise) 

t section pitching (or torsion) moment per unit length 

perpendicular to plane of symmetry 

R concentrated pitching moment or torque 






c tan A 


(EI)j. cos a 
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S 

Vt 

y 

y* 

a 

a 

Aa 

a© 

r 

8 


6 

11 * 



A 

\ 

«P 


wing area 
volume of tip tank 

lateral ordinate measured from plane of symmetry 
dim ensionless lateral ordinate 


local angle of attack, radians 

total angle of attack, including increment due to aeroelastic 
action, radians 

angular deformation of iDOom at lifting surface, radians 
effective angle of attack due to unit aileron deflection 


/dC2/dS \ 
\^cijdal 


local dihedral or spanwise slope of elastic axis, radians 

aileron deflection (in planes parallel to plane of symmetry), 
radians 

moment-arm ratio (®2/®l) 
lateral distance from wing root 


dimensionless lateral distance 



angle of sweephack at elastic axis 
taper ratio (Tip chord/Root chord) 
angle of twist about elastic axis, radians 


Subscripts: 

a pertaining to antisymmetric case or to aileron 


ED 


at divergence of boom 
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D at divergence 

g geometric (built in or due to airplane attitude) 

P due to concentrated normal force 

R at reversal or due to concentrated pitching moment 

r at wing root or reference value 

s due to structiiral deformation 

t at wing tip or pertaining to tip tank 

V pertaining to boom-mounted lifting surface, or in the 

presence of the lifting surface 

a due to angle of attack 

8 due to aileron deflection 

A referred to axes parallel and perpendicular to elastic axis 

0 rigid wing = "q = 0) 

Superscripts: 

c due to concentrated force and moment 

P diie to concentrated normal force 

R due to concentrated pitching moment 

8 due to aileron deflection 


Matrices: 


[ ] 

square matrix 

{} 

column matrix 

L J 

row matrix 

L 1 

diagonal matrix 
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[B] 


idjentity matrix 

matrix defined "by equation ( 29 ) 
aeroelastlc matrix (equation (l3a)) 
auxiliary aeroelastlc matrix (equation (l3t>)) 
aileron-reversal matrix (equations (2l) and ( 32 )) 
matrix defined by equation ( 28 ) 


{h}, {h} 


matrices defined "by equations (l4) and ( 30 ) 


[l],[ll] single and dOTible integrating matrices from tip to root 

(prime mark on symObol I or II designates an Integrating 
matrix for a function which goes to zero with Infinite 
slope at wing tip) 


f I 

[I] single integrating matrix from root to tip 

[l^J, f’irst rows of matrices ^1] and [II] 


[Q] matrix, of aerodynamic Influence coefficients 


METHOD OF ANALYSIS 
Matrix - Integration Method 


Resume^ of method of references 1 and 2 . - The method of references 1 
and 2 Is based on numerical integrations of the equations of structural 
equilibrium by means of suitable integrating matrices. These Integrating 
matrices, together with other matrices, and constants which describe the 
structural, aerodynamic, eind geometric properties of the given wing, are 
assembled into aeroelastlc, auxiliary aeroelastlc, and aileron-reversal 
matrices, from which the structural twist at any dynamic pressure, as 
well as the dynamic pressures required for divergence and reversal, can 
be determined. The method of references 1 and 2 was modified slightly 
in the following resume. 

The limitations of the method of references 1 and 2 are discussed 
in those papers; they arise from the assumptions that the spanwise lift 
and pitching-moment distributions csn be predicted for any given angle- 
of-attack distribution by means of the aerodynamic influence coeffi- 
cients and that the struct\iral deformations can be predicted by simple 
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beam theory plus rigid-body rotations imparted by the root, (in refer- 
ences 1 and 2 a method Is also presented for using structural Influence 
coefficients measured on the actual wing or calculated by methods more 
refined than simple beam theory; this method can be extended to the 
case of wings with concentrated forces at the tip in the same manner as 
employed in the present paper for the method based on sin 5 >le beam 
theory. ) 

The loading coefficient cc^/'c for any section of the wing 

parallel to the stream may be determined for any angle- of-at tack dis- 
tribution by means of suitable aerodynamic Influence coefficients Qg 

and Qa (for symmetric and antisymmetric lift-distributions, respec- 
tively) in the form 

L. & 


where a is the total angle of attack at a given point on the span due 


to geometrical setting and structural deformation, 
wing lift-curve slope. 


"La 


is the rlgid- 


Cj is the negative of the coefficient of 
‘'d 


damping in roll, and •< 


CCi 


cCi 


is 


times the section loading coef- 


J 5 


ficient due to aileron deflection for a unit equivalent angle of 


attack ctg6; the matrix 


Approximate Influence- coefficient matrices Q 

lated for subsonic flow by the method of reference 3* 
section can then be written as 


is used for the ssike of definiteness. 

and may be calcu- 

The lift on any 


{1} - 1C [S|£j. (2) 

and the section pitching moment about the elastic axis in planes 
parallel to the free stream can be written as 






( 3 ) 


where the subscripts a and 6 serve to specify the lifts due to the 
loading coefficients represented by the first and second term on the 
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ri^t-hand side of equation (l), respectively, and wliere 6 is the 
ratio of the moment arm B 2 to the moment arm e^L. (See fig. 1.) 

The parameter eQ_^ is an arbitrary reference value of the dimension- 
less section moment arm e-]_. 

The section lift and section torque given by equations (2) and (3) 
can be Integrated by means of Integrating matrices to obtain the 
accumulated bending and twisting moments at any section about a pair of 
axes parallel and perpendicular to the free stream, respectively, with 
their origin at the elastic axis at that section. Thus, 



(4) 

= ^DO {t> - tan A {m} 

(5) 


where the matrices [Y] and are defined and given in refer- 

ence 1. If the lift distribution goes to zero with infinite slope at 
the wing tip, the matrices must be modified to take this fact into 
account; the resulting matrices are designated by [^*3] and 
respectively. 

The bending and twisting moments obtained in this manner can be 
transferred to axes along and peipendicular to the elastic axis; the 
resulting moments are 


= cos A - sin A 

(6) 

^T^ = sin A -{mJ- + cos A -[tJ- 

(7) 


The structural twist cp and the slope of the structural-deformation 
curve r can then be obtained by an Integration of the products Ty^/GJ 

and M/JeI, respectively. These integrations can be performed numeri- 
cally by means of the matrix ’s^iich, for equally spaced stations, 

is the double transpose of the matrix [Y]. Consequently, 



^1/^ fil't 

(GJ)j. cos a* — 



( 8 ) 
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{hj. - [a] 



(14) 


^*/2 tan A 
(EI)^ ®1j.^ cos^A 


( 15 ) 



( 16 ) 


and where Oq Is the effective angle of attack due to unit aileron’ 
deflection. If the rlgld-wing, lift-curve slope Cj^ in equation (l2) 
is based on a nonlinear lift curve, the value of Cj^ should be taken 

at an average angle- of-attack condition. For symmetric cases the 
matrix is used Instead of in equation (13) and the second 

term on the ri^t-hand side of equation (ll) is disregarded; for lift 
distributions ■which go to zero with infinite slope at the tip, the 

matrices [l*]] and [jl*J 3-re used in equations ( 13 ) and (1^) instead 

of [l] and [i:Q . 


The aeroelastic equation (ll) can be solved for cCg and for the 

dynamic pressure at divergence or reversal in the manner described in 
references 1 and 2. Specifically, the reversal speed can be obtained 
by determining the eigenvalues of a matrix obtained in the fol- 

lo’wlng manner; At reversal the rolling moment due to aileron deflec- 
tion is to be zero, so that 






( 17 ) 


where 




( 18 ) 
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and where, in turn, 
matrices [ll^ and 


[Ilij a»d [Iij 
[l] . Hence, 


are the first rows of the 


Og8 


L“oJ [aj 




CCj 

< _ 
cCt 
^ d 

to. — • 

8 


( 19 ) 


so that, upon suhstitutlon of eq.uatlon ( 19 ) into the second term on 
the ri^t-hand side of equation (ll), this eq.uation becomes 


-[CC3J. 

where the aileron- reversal Diatrix JAjQ is defined by 

[Ah] = M + I W M 

and 



( 20 ) 


( 21 ) 


( 22 ) 


In the derivation of eq.uation (20), the total angle of attack has 

been replaced by in equation (ll) because the geometric angle 

of attack has no effect on aileron reversal. 

Modifications required for inclusion of tin forces .- As a result 
of the fact that the method of references 1 arid 2 uses the equation of 
structural equilibrium in Integral rather than in differential forma 
the inclusion of concentrated tip forces is accomplished quite readily 
by including additional matrices which introduce the effect of the con- 
centrated forces at the tip in the aeroelastlc and aileron-reversal 
matrix. These matrices are then treated in the same manner as those 
for the wing without concentrated forces. In essence this procedure 
amounts to performing a separate analysis for the wing with and without 
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■tlie concentra'bed forces, alttLOU^ many of iJie matrices calculated in 
the wing-alone analysis can be used in the other analysis. This method 
is subject to the same limitations as the method of references 1 and 2 ; 
in particular, the aerodynamic forces on the wing must be predictable 
by means of suitable Influence coefficients. 

If the normal force at the tip is and has a pitching moment T^ 

about the elastic axis in a plane parallel to the plane of symmetry, 
equations ( 4 ) and ( 5 ) become 


CM} - 

(23) 

= ^[l]-Ct> - tan A-CM> + {t^J 

( 24 ) 


The column matrices 



consist of elements all equal 


to and T^, respectively, and Is the dimensionless lateral 

distance from the root to the station at which the bending and twisting 
moments are obtained. If the concentrated force and moment are due to 
aerodynamic action they may conveniently be expressed in terms of dimen- 
sionless coefficients as 


or 


^t = 

Tt = CLo^q,dS^((X(. + KS) 


T4. = 


(25) 

(26a) 

(26b) 


where oL|. Is the value of the total angle of attack at the tip and 
where equation (26a) pertains to a boom-mounted lifting surface with a 
center of pressure at a distance d ahead of the el^tic axis, and 
equation (26b) pertains to a tip tank with volume V^. For the sake of 

definiteness equation (26a) is used in the following analysis. The 
factor K in equations (25)# (26a), and (26b) is the gear ratio between 
the boom deflection and the allercm deflection. In the case of a tip 
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tank or In the case of a lifting surface which is not geared to the 
aileron, the factor K is zero. 

When equations (4) and (5) are replaced by equations (23), (24), 
(25)^ and ( 26 a), the aeroelastic equation (11) becomes 

{ag} = q*|[A] + [:B][]^.]Jfa> - q.%5 {h} (2?) 

where 


[B] = 




S+ 


c(bV2) 


[I] 


t t 


d 

|(GJ)r] 

(GJ)r 

■^r 

Lgj 

(El)r 


P 

tan A 





( 28 ) 



0 0 0 1 
0 0 0 1 
0 0 0 1 
0 0 0 1 


(29) 


W ■ {»} - ^[b]© 


(30) 


and -{]l^ is a column all the elements of which are equal to 1. 

The aeroelastic equation (27) can be solved in the same manner as 
equation (11) for or for the dynamic pressure at divergence. An 

aileron- reversal matrix can be calculated in the manner employed to 
obtain equation (22). The rolling moment due to aileron deflection 
vanishes when 
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Ct 

_St 

c(V2) 


KB + 


St 

°Lq_ c(b/2) 



0 


( 31 ) 


where the last two terms on the left-hand side represent the rolling 
moment of the lift on the boom-mounted lifting surface and where the 

matrix |_ltj first row of the matrix . If otgB is 

obtained from equation ( 31 ) substituted into equation (27) the 

resulting equation is identical in form with equation (21), except that 

the aileron- reversal matrix is replaced by a new matrix 



-DQ + H[h] -^{5} 


\Foj [«a] 


+ 



St 

c(b/2) 



(32) 


where 


g’ 





Ct 

cti. Si xr 


( 33 ) 


Simplified Method 

The simplified method is applicable in cases for which the aero- 
dynamic interaction between the wing proper and the source of the con- 
centrated aerodynamic force at the wing tip is neglected. This method 
consists of determining the aeroelastlc twist of a wing subjected to 
concentrated torques (or pitching moments) and normal forces of known 
magnitude applied at its tip at one dynamic pressure. This information 
is combined with the aerodynamic characteristics of a boom-mounted 
lifting surface, and the results are extrapolated over the range of 
dynamic pressures of interest on the basis of the semi rigid- wing con- 
cept; the dynamic pressures at divergence and reversal are determined 
from this extrapolation. 
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This sinrplifled method is subject to the same limitations and to 
two others, as well. The semirigid concept furnishes a useful basis for 
extrapolation of aeroelastic results only if the lowest root of the poly- 
nomial for the dynamic pressure at divergence is much lower in absolute 
value them the next hi^er one, as seems to be the case for actual wings. 
Also, as developed in this paper, aerodynamic interaction between the 
source of the concentrated force and the wing proper is not taken into 
account. If the magnitude of this interaction can be predicted it can be 
taken into account in the sinplified method by certain modifications, as 
discussed in a later section, but in such a case it may be more e3q>edient 
to use the other method presented in this paper. 

The effect of concentrated forces of known magnitudes on the aero- 
elastic characteristics of a wing alone .- If a concentrated normal 
force P and a concentrated pitching moment R are applied to the wing 
tip, the bending and twisting moments about axes parallel euid perpen- 
dicular to the plane of symmetry are 


{m} = p ^ {1 - (34) 

and 

^Tj- = -tan A + R -£l}’ (35) 


From these moments the angle-of-attack change Og*^ caused by these con- 
centrated forces can then be calculated from equation (10) by using 
equations (6), (7), (8), and (9). The resulting expression for OgC 
may be written as 



R 


cos A(b«/2) [-,-1 
(GJ)r ^ 


t I 


(GJ)r 

GJ 



+ 


sin^A(b*/2) r^-i 1 1 
(EI)^ cos a 



P 


tan A(b*/2)* 
(El)y cos a 


■M 


I I 




(36) 
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where is the angle- of-attack change produced directly by the con- 

centrated forces without the presence of aeroelastic effects. If the 
concentrated forces are caused by a boom-mounted lifting surface and if 
the aerodynamic- Induction effect of the wing on the lifting surface can 
be neglected, then 


and 


^ ( 37 ) 

E = CL^S^q.d(o 4 . + K5) ( 38 ) 


where 04 . is the angle of attack of the wing at the tip and Includes 

both the angle of attack due to airplane attitude and that due to 
structural deformation. Substitution of equations (37) and ( 38 ) into 
equation ( 36 ) yields an escpression which may be reduced to 


= ('^ + K5) 1* 



(39) 


where 


Kt> = 


■— Cty' 


Kr = 


a 

— °v d 






(40) 


(41) 


( 42 ) 
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{a®}. [I]” 


(GJ)r 

_GJ 


+ 


(GJ)r 

(EI)^ 


p 

tan A 


(El)r' 


El 





m 


or, if 


GJ _ El 
(GJ)r ~ (El)r 


equation (^3) reduces to 




(GJ)r 

(EI)^ 


tan^Al 



m 


Aeroelastic effects on the angle- of-at tack change due to concen- 
trated forces .- The aeroelastic effects of the angle- of-attack change 
due to the concentrated forces given by equation (39) can be calculated 
by introducing this angle-of-attack change in the rl^t-hand side of 
equation (11) provided that the aerodynamic- Induction effects of the 
boom-mounted lifting surface on the wing are neglected. Hence 


{°^s} = J {agj. + {ae} 




{h} 


or 


U1 - q.*[Aj {o-aj = q*[A] + q*[A] - q*aQS ^h}- (45) 


where column which describes the angle-of-attack changes 

caused by aeroelastic action due to all three forcing loadings or angles 

of attack {»g). and •^hj’ aond ■vAiere, in turn, ts'} can be 

considered to consist of two parte, as indicated in equation (39)* The 
most convenient way of solving equation (45) consists of evaluating 

separately the contributions of this purpose 

equation (45) can be rewritten as 
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[111 - <1*[a1] {Sgj- - {ag} 

(h6) 

[Ul - 4*[Al] M - {=‘^} 

ihl) 

[Ul - 4*[Al] {a®} = {a?} 

(48) 

[Ill - a»CA]J " "“s® ^**3- 

(i^9) 


wliere "a_ is ecjual "to a_ plus tlie part of a_ due to a_ and thus 

is the total (or net) angle of attack due to airplane attitude and the 
amount of aeroelastlc deformation associated with that angle of attack 

at the given value of <i*; similarly ^ and' c? are eq.ua! to a? 
and oP plus the amount of aeroelastlc deformation associated with 
these angle-of-attack distributions at the given value of q.*, and 
when multiplied 'by q.*, is the amount of aeroelastlc deformation associ- 
ated with aileron deflection. In the matrix [a] -(dilch occurs in eq.ua- 
tlons ihS), and (48), [Qg] has to be used for symmetrical flight 

and for anti symmetrical fll^t; In eq.uatlon (49), is used. 

The total ■ angle-of-attack distributions due to all forcing angle-of- 
attack distributions and their associated aeroelastlc increments are then 

{a> = jkg} + (at + K6)q* jKR|a^Jy - Kp + q* ^ag^J. (50) 

and hence, at the wing tip 


^^^t = “&t (“t ^ - KpOt^) + q*Ug 


(51) 
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so that the angle of attack of the bocan-mounted lifting surface is 


+ KS = 


K 6 4- 

1 - - KpOt^) 


(52) 


The static aeroelastic analysis of a ving with a boom-mounted 
lifting surface at a given value of q.* may therefore be performed as 

follows: The columns obtained as part of the 

analysis of the wing alone. For the vane, the parameters and Kj^ 

as well as the columns {g,P} and {g,®} are calculated from equa- 
tions (4o), (^l), (42), and (43). Hence, the columns 
are obtained by solving equations (47) and (48), respectively. ^If 
Grout’s method (reference 4) was used to obtain ®uid then 

the evaluation of the two new columns requires very little additional 
effort; if series-expansion or iteration method was used to calcu- 
late ^ then a new iteration is required, which should 

converge more rapidly than that for 'j^^gj*' > sii^ce experience indicates 

that and tend to approximate the dominant modal column 

of the matrix [XI more closely than does "^g^ angle of attack 

of the lifting surface can then be determined from equation ( 52 ), hence, 
the new angle- of-attack distribution from equation ( 50 ) and, finally, the 
lift distribution from equation (l). The lift on the lifting surface 
itself can be obtained from equation ( 37 ) • 


Extrapolation of aeroelastic characteristics calculated for one 
value of q* to other values of q* .- The foregoing procedure can be 
repeated for any value of q* of interest, but the values of and 

q* cannot be obtained directly from this analysis. In order to calcu- 

X\ 

late these values and to permit the extrapolation of results calculated 
by the method indicated in the preceding paragraph, the semirigid-wing 
concept may be used provided that the lowest root of the polynomial for 
the dynamic pressure at divergence is much lower in absolute value than 
the next hi^er one. In essence, this concept consists of reducing the 
degrees of freedom inherent in the structural deformations of the wing 


KAGA EM L52A22 


21 


to two "by stipulating shape of the bending and twisting deformations 
and calculating the magnitude of eadi. 

For the purpose at hand, the results of a semirigid analysis can 
he obtained by considering a rigid constant- chord wing permitted to 
rotate about hinges at its root parallel and perpendicular to its 
leading edge subject to the restraint of a torsion and a bending spring 
with constants and respectively. In this case, the lift on 

one half-wing is 


I = l(“g + M 

The twisting moment about the torsion hinge is 

T^ = e^c cos A ^ 


The bending moment about the bending hinge is 



h/2 L 
cos A 2 


where y* is the dimensionless lateral center of pressure. The angle 
of twisting deformation is 


cp - K^T^ 


and the angle of bending deformation is 




so that the angle of attack due to structural deformation is 
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ttg = cp COS A - r sin A 

= q*(l - k)(ag + ttg) (53) 


where the parameters q* and k are similar to those previously 
defined but are defined in equation (53) as 


q* = qKrj cos^APj^ej^c | 


J5. ^ _M tan A y*(h/ 2 ) 
%' cos^A ®1^ 


The solution of equation (53) can be expressed either as 


°^s = 


< 1*(1 - k ) 

1 - (1 - k)q* 


ct, 


■g 


or 


a = 


1 - (1 - k)q» 


cx, 


g 


where a = ttg + ag. Since divergence will occur when 


= 


1 - k 


the preceding equations can be written as 


oc„ = 
B 


q*/q^ 




a 


g 


( 5 ^f) 


1 
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and 


a 



( 55 ) 


As shown In reference equation (5^) yields a good approximation 
to the angle of structural deformation of an actual wing for all values 
of q*, if a constant C obtained from an analysis for one value of q*"^ 
is introduced. This constant is different for each point along the span 
and for each geometric angle- of-attack condition. With this modification 


a 


= C 


q*/q*jj 



a, 


■s 


( 56 ) 


where q* and q*jj are now defined in accordance with equation (12) 
for the given wing. Hencej 


a = 


1 - (1 - c)(q*/q.*p) 


1 - 


Q.* 

<l\ 




( 57 ) 


For the wing with a boom-mounted lifting surface the following 
approximate expressions can be written: 


1 - (1 - 


CL^ — 

St 


1 - 




a, 


®t 


( 58 ) 


St^ = 


1 - (1 - C2)(qV<l*D) _ R 

Ctj. 


1 - 


q.* 


( 59 ) 
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p ^ 1 - (1 - C3)(q.M) p 

1 < 1 * 


( 60 ) 



6 = -agS 


1 


- (1 - 





( 61 ) 


■wliere q*jj is the value at divergence of the parameter q* defined by 
equation (12) for the wing alone and -where hj. is the last element of 

the column ■^h^. By calculating , and at one 

value of q* from equations (46), (47), (48), and (49), the constants 
C 3 _, Cg, 0^, and Ci,. in equations (58), (59)^ ( 60 ), and (6l) can be 

evaluated. These equations can then be substituted into equations (52) 
and ( 50 ) to yield, respectively. 



Equation (62) gives the value of c^. + K& at any value of q-*^. 
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Calculation of q* at divergence and reversal .- The value of q* 
at divergence of the wing with 'a boom-mounted lifting surface can be 
obtained by equating the denominator of the fraction on the rl^t-hand 
side of equation (62) to zero. This procedure yields the quadratic 
equation 






- - (1 - C3)KpC4.P) = 0 


(64) 


which can be solved for 



and, hence, for 



the value of 


q* 


at divergence of the wing with a lifting surface. Of the two value 
of q* obtained in this manner, the smaller one in absolute magnitude 
■^v 

is the critical one unless It corresponds to a negative value of 

in which case the larger is the critical one unless it also corresponds 
to a negative c^. In which case the wing cannot diverge. 


In order to calculate the reversal speed of the wlng-wlth- lifting- 
surface combination, the sum of the rolling moments of the lift distri- 
bution due to structural twist and to aileron deflection and of the 
rolling moment caused by the lift on the boont-mounted lifting surface 
is set equal to zero, as in equation ( 31 ) which may be rewritten as 
follows ; 


L^^oJ [<^]{“> + 



+ + K5^ 


= 0 


(65) 


where 


«v = 


c(b/2) 


( 66 ) 
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The value of required In equation ( 65 ) can be obtained from equa- 

tion (63). However, the effect of the approximations made in obtaining 
equation (63) can be minimized by first substituting equation (50) into 
equation (65) and then approximating the moments o_f the various lift 
distributions in the manner employed to approximate the angles of attack 

in equations (58) to (6I). For the sake of convenience -^aj- may be 

calculated for ccg6 = 1 and OgS set equal to 1 in equation (65); in 

the subsequent derivation this simplification is assumed to have been 
made. This procedure yields the equations 


L^IoJ M |kr {a?} - Kp + q* 


> + 






or 


{•H * - (1 - ‘=’ 2 )^) - * 


- $;)) - ^ - %) 


= 0 


(67) 


where Kq, Kp, K^, K^,., C’p, *^'3# defined by the 

relations 


L^^oJ 


CCi 


cC 


^d 


> =Ko 

5 


( 68 ) 
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1 - 

L^ioJ M ^ 3* - 

■ i*D 


L“oJ M (4 




( 69 ) 


(TO) 


( 71 ) 


The coefficients K2^ K^, and Kij. are eq,ual to the matrix products on 

the left-hand sides of equations (69) j (70), and (71) with the 

columns £h^ substituted for 1^°^' 

anil respectively. The coefficients 0*2^ 

then he obtained by evaluating the left-hand sides of eq,uations (69), 
(70), and (71) with the columns obtained by solving equations (^ 7 ), (^), 
and (49) at one value of <1* and substituting that same value of 4* 
and the previously calculated values of Kg, K^, and K|j^ on the rl^t- 

hand sides of equations (69), ( 70 ), and ( 71 )- 

The value of (X|. + Ky/og given by equation ( 62 ) can be substituted 
into equation (67) to yield 
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By collecting the terms of equal powers of , equation (72) can be 

reduced to a fourth-degree polynominad in q* . Of the four roots of 

this polynominal the lowest real one of appropriate sign (to correspond 
to a positive value of q) represents the critical aileron reversal 
speed. 

Calculation of the effect of boom flexibility .- The effect of the 
bending flexibility of the boom on static aeroelastic phonomena can 
easily be taken into account in the method of the preceding sections. 

The flexibility of the boom may be defined by the angle-of-attack change 
of the lifting surface due to boom deformation per unit normal load 
applied at the lifting surface in still air Kg. 


The change in angle of attack of the lifting surface due to boom 
flexibility is then, as a result of equation ( 37 )^ 


or 


Zi&t = EKg 


= + kb + 


(73) 


04 . + K5 + + KS) 


1 - 


(74) 


q*. 


BD 


where Q.*gj) "the value of q* required for divergence of the lifting 

surface as a result of boom flexibility, tliat is, considering the wing 
rigid. The value of is given by 


d ^ 


cos A 




Kb^r(GJ)j 


(75) 


Equation (74) indicates that in order to take boom flexibility into 
account the angle of attack of the lifting surface 04 + K5 must be 
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replaced by the product of this angle and the factor 



This 


procedure is equivalent to multiplying either the lift-curve slope or 
the area of the lifting surface by this factor. Consequently^ equa- 
tions (50) i C52), (62), and (63) are valid for the flexible-boom case 
provided the factors Kp and Kp in these equations are divided by the 
a* 

factor 1 - — = — . The angle of attack of the lifting surface relative 

to the free stream can be obtained by dividing the values given by equa- 
tions (52) and (62) by this factor. The dynamic pressure at divergence 
can be obtained from equation ( 64 ) if the term l/q*^ is added to the 


three terms -vd-thin the second parentheses and to the two terms within 
the third parentheses on the left-hand side of that equation. Similarly, 
the dynamic pressure at reversal can be obtained from equation (72) if 


the term 





underlined in equation 


the factor 




( 72 ) 


is multiplied ly 


The procedure outlined in this section can also be used to take into 
account the effect on static aeroelastlc phenomena of the flexibility of 
the lifting surface Itself by calculating the value of q* required to 
diverge the lifting surface and using this value instead of 


Comparison of the Two Methods of Computation 


A conqjarison of the numerical results of the two methods may be had 
frcm the following tabulation of some of the results for the case of the 
wing - lifting- surface comblnatian at subsonic speeds with = 0 . 02 , 

— = 1 . 5 ^ gear ratio 1 , and an infinitely stiff boom discussed in a 

subsequent section of this paper: 


Method 


q*j^( 3 ) 

Matrix- integration 
Simplified 

-0.1692 

-.1699 

1.368 

1.375 
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In this tahalation q* (3) is the third root (in absolute magnitude) of 

the polynomial for "the lowest two roots are conqplex conjugate 

numbers and, hence, have no physical significance. There is good agree- 
ment between the results of the two methods. 

In view of the satisfactory agreement of the results of the two 
methods and in view of the fact that the simplified method is generally 
less time consuming than the matrix- integration method, the simplified 
method appears to be preferable in all cases where it is applicable, 
particularly when wing-alone calculations have been made previously or 
when a number of configurations Involving different tip forces are to 
be analyzed for the same basic wing. The matrix- integration method is 
more widely applicable than the simplified method; when the simplified 
method is applicable the matrix- Integration method is preferable only 
in the case Tdiere the source of the concentrated load is permanently 
installed, so that no wing-alone calculations need be made. 

The extension of the two methods presented herein to the calcu- 
lation of the aeroelastin effects of concentrated forces located at 
points on the span other than the tip presents a problem in that such 
a force gives rise to moment and torque distributions which are either 
discontinuous or have a discontinuous slope. Such distributions cannot 
generally be integrated accurately by the siiip>le numerical methods on 
which the Integrating matrices used in this paper are based. However, 
special integrating matrices which take these discontinuities into 
account can be set up for the purpose of calculating the structural 
deformation for any concentrated force or moment at a given point on 
the span. Also, Interpolating matrices can be devised for calculating 
the angle of attack at a given point on the span in terms of the angle 
of attack at the points on the span used in the aeroelastlc analysis. 

By incorporating these interpolating and special Integrating matrices 
in the method of references 1 and 2 in a manner similar to that indi- 
cated for tip forces in the present paper, a method can be obtained for 
taking concentrated aerodynamic forces at points other than the tip 
into account in aeroelastlc calculations. 


SCOPE OF THE ILLUSTRATIVE CALCULATIONS 
Unswept Wing with Tip Tank 


The matrix- integration method presented in the preceding section 
has been used to calculate some static aeroelastlc characteristics of 
an unswept wing with a tip tank. The geometric and some of the struc- 
tural and aerodynamic characteristics pertinent to the aeroelastlc 
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analysis as well as a plan, form of tlie tip tank are given in table l(a). 
The wing plan form and, the tank are the same as those used in refer- 
ence 6, and some of the aerodynamic data used in the present paper have 
been obtained from that reference; the wing lift and moment coefficients 
are those for the wing with section B with and without the tank with 
sealed gap at a Mach number of 0.8. The tank lift and moment coefficient 
(referred to the mean aerodynamic chord and one-half the wing area) are 
for the tank on the wing with section A and gap open, since these data 
are not available for section B and gap closed. The lift-curve slope 
and the moment-cuive slope given in table l(a) for the taink-on configu- 
ration are those for the wing alone in the presence of the tank; they 
have been obtained by subtracting the lift and moment on the tank from 
the total lift and moment on the wing-tank combination. 

The stiffness distributions El/(El)j. and Gj/(CrJ)j. of this wing 

are assumed to be Identical and are given by the dashed- line curve in 
figure 2. They were obtained by means of the constant-stress concept 
of reference 5 for the inner 70 percent of the semi span; in the outboard 
30 percent of the semispan they are assumed to vary as the fourth power 
of the chord. 

The rigid-wing spanwise lift distributions of the wing alone, for 
uni form angle of attack, for linear antisymmetric twist and due to 
aileron deflection were obtained by the method of reference 7 an<i are 
given in figure 3 tiy the lines labeled q.* = 0. The spanwise lift 
distribution for the wing with the tank at a unlfom. angle of attack was 
estimated by distributing the additional lift carried by the wing due to 
the presence of the tank near the tip; the resulting distribution over 
the part of the wing not blanketed by the tank is shown in figure 3Ca-)- 
The rigid-wing lift distributions for the other two angle- of-attack con- 
ditions were then estimated by using the method of reference 3 fn. con- 
junction with the lift distribution for uniform angle of attack estimated 
in this manner. The factors k]_, kg, k^, and k|^ required in the 

method of reference 3 were obtained from the figures of reference 3 for 
the aspect ratio which a wing without a tank would have to have in order 
to have the same lift- curve slope as the actual wing in the presence of 
the tank. The rigid-wing lift distributions for the part of the wing 
not covered by the tank calculated in this manner are shown in fig- 
ures 3(b) and 3(c) by the lines labeled q* = 0. Aerodynamic Influence 
coefficients for this wing were calculated by the method of reference 3 
using the lift distributions shown in figures 3(a-) and 3(t)). 

The spanwise variation of the local aerodynamic- center positions 
of the wing alone was estimated from an analysis of lifting-surface 
calculations and experimentally obtained pressure distributions on 
similar wings and was adjusted to correspond to the pitching moment 
measured in reference 6. This variation was modified sll^tly for the 



32 


MCA m L52A22 


tank-on configuration by using the assumed spanwise variation of the 
additional lift carried by the wing due to the presence of the tank as 
well as the assumed chordwise location of this Increment in lift obtained 
on the basis of the pitching-moment data of reference 6. The spanwise 
variation of the moment arm e^, which is the difference between the 

local aerodynamic- center and elastic-axis locations, is shown in fig- 
ure 4(a) . The reference value of e^^^ was taken as the value 

of wing alone. 

The local centers of pressure due to aileron deflection were 
obtained from the assumed section center of pressure due to aileron 
deflection (42 percent chord), the local aerodynamic- center positions, 
and the spanwise lift distributions due to aileron deflection by the 
method outlined in reference 2. The dimensionless distances e 2 of the 

centers of pressure due to aileron deflection from the elastic axis (see 
(fig. 1) are also given in figure 4(a). 

In the aeroelastic calculations the wing was assumed to be mounted 
on a reflection plate, as in the tests of referenice 6, so that b* = b. 
The small angle of sweepforward of the elastic axis (1.6°) was neglected. 
All root-rotation constants (see reference 1) were assumed to be zero. 

Calculated for wing with tank on and off were the dynamic pressure 
at divergence and the dynamic pressure at reversal; also calculated for 
several values of were the spanwise lift distributions due to 

uniform angle of attack, due to linear antisymmetric twist, and due to 
aileron deflections, the lift-curve slopes, the coefficients of danping 
in roll, the rolling-mcment coefficients due to aileron deflection, the 
spanwise centers of pressure, and the rolling velocity per unit aileron 
deflection- 


Sweptback Wing with Boom^Mounted Lifting Surface 

The geometric characteristics, as well as sane of the assumed aero- 
dynamic and structural characteristics of a 45° sweptback wing, for 
which aeroelastic calculations similar to those described in the pre- 
ceding section have been performed, are presented in table 1(b). The 
stiffness distributions have been estimated in the same manner as that 
employed for the unswept wing and are presented in figure 2. The rigid- 
wing spanwise lift distributions at subsonic speeds were calculated by 
the method of reference 7 and are shown in figure 5 i^he lines 
labeled q* = 0. Aerodynamic Influence coefficients for subsonic speeds 
were calculated by the method of reference 3. For supersonic speeds 
strip theory was used; the resulting lift distributions are shown in 
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figure 6. Tlie moment arms e^_ and for sutsonlc speeds wBre esti- 

mated in the manner employed for the unsvept wing and are shown in fig- 
ure ij-Ch); for supersonic speeds the values were estimated from linearized 
two-dimensional theory. Reference values for e-|_^ of 0. 2 and 0 were 

used arbitrarily in the calculations for subsonic and supersonic speeds, 
respectively. 

As in the case of the unswept wing the sweptback wing is considered 
to be mounted on a reflection plate, and all root- rotation constants are 
assumed to be zero. The aerodynamic interaction between wing and boom- 
mounted lifting surface has been neglected in the calculations. 

No specific boom-mounted- lifting-surface plan forms have been con- 
sidered; the surfaces are characterized in the calculations by the ajrea 
ratio 




by the moment-arm ratio by fbe gear ratio K of lifting- surface 

motion to aileron motion, and by the boom flexibility Kg or the dimen- 
sionless dynamic pressure for boom divergence <l*gg defined in equa- 
tion (75) • Calculations have been made for the combinations of these 
parameters shown in table 2. The combinations for which ^ = 0 have 


no physical significance and are used only to Illustrate certain trends. 


For each of these combinations the aeroelastlc Information listed 
at the end of the preceding section was calculated using the simplified 
method; for most of the combinations, excluding those with flexible 
booms, calculations were also made by the matrix -integration method. 

For a configuration with = 0.02 and — = 1.5 calculations have 

been made for an ungeared lifting s^lrface (case 3) as well as for a 
geared lifting surface with gear ratio K = 1 (case 4); the calcula- 
tions for all other cases have been made only for a gear ratio of 1. 
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RESULTS AND DISCUSSION 
UnBvept Wing with Tip Tank 


The. lift distrihutionB of the unswept wing with and without tip 
tank for the three angle-of-attack conditions considered are shown in 
figure 3 at sevei«.l values of q.*. The values of q* for the wing with 
the tank are based on the value of Cj;^ for the wing without the tank, 

so that for the same value of q* the dynamic pressures for the wing 
with and without the tip tank are the same, since all other quantities 
that enter into the definition of q* are the same for both cases. 

The effect of aeroelastic action is, as eaqiected, to increase the 
lift at all points on the span, particularly in the region near the tip. 
This increase is much more pronounced for the wing with tip tank than 
for the wing withxjut tip tank; even at somewhat lower values of the 
dynamic pressures (q* = O.I 92 as opposed to q* = 0.255) the increase 
in lift on the wing with the tip tank is much greater than that on the 
wing without the tip teuik. These two values of q* represent the same 
fraction of q*j^ and differ from each other because q*j^ is different 

for the two cases. 

The wing lift coefficient, lateral center of pressure, rolling- 
moment coefficient, and rate of roll obtained by integrating the lift 
distributions shown in figure 3 are presented in figure 7 as functions 
of the dimensionless dynamic pressure q* preferred to the value of 

for the wing without tip tank, as in figure 3 ^- lift is seen to 

Increase much more rapidly for the wing with the tip tank than for the 
wing without the tip tank; the spanwlse center of pressure is farther 
outboard at q* = 0 and moves outboard more rapidly with increasing q* 
for the wing with the tip tank than for the wing without the tip tank. 

At q* = 0 (rigid wing) the rolling-moment coefficient due to unit 
aileron deflection is 0. 220 for the wing without tip tank and 0 . 29 I for 
the wing with tip tank; for the wing without tip tank it decreases with 
increasing q*, whereas for the wing with tip tank it Increases with 
Increasing q*. The coefficient of daji^tlng in roll is 0.436 for the 
wing without tip tank and O .685 for the wing with tip tank; it increases 
with increasing q* in both cases but much more rapidly in the case of 
the wing with tip tank. The rate of roll is less at q* = 0 for the 
wing with tip tank than for the wing without tip tank and decreases more 
rapidly with Increasing q*. 

The value of q* required to diverge the wing without tip tank is 
1.021, and the value for the wing with tip tank is O. 38 O; the value 
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of q,* required to reverse the lateral control of the wing without tip 
tank is 0 . 819 ^ whereas for the wing with tip tank the reversal speed is 
higher than the divergence speed, the value of q*jj being 0.409. (The 

value of q* for antlsynunetric divergence of the wing with tip tank 

is 0 . 388 .) 

As may be seen from the results presented in the preceding para- 
graphs, the tip tank tends to have a very unfavorable effect on the 
static aeroelastlc characteristics of the wing; for Instance, the 
dynajnic pressure required to diverge the unswept wing with the tip tank 
considered in this paper is very much lower than that required to 
diverge the wing alone. This IB due in part to the higher lift carried 
by the wing and the more forward local aerodynamic centers, particularly 
near the tip, that result from the presence of the tank and in part to 
the concentrated moment introduced by the tank proper. Consequently, a 
wing which does not diverge by itself may diverge in the presence of the 
tip tank. Actually, the wing may destroy Itself even before reaching 
the lowered divergence speed, because as it approaches this speed, the 
lateral center of pressure moves so far outboard and the lift-curve 
slope becomes so large that a relatively small gust may overstress the 
wing. 


The values given here for the decrease in dynamic pressure required 
for divergence and increase in severity of the static aeroelastlc 
phenomena, in general, may be somewhat pessimistic for two reasons. The 
stiffness distribution assumed for the wing is likely to be too low near 
the tip compared to actual airplanes. A somewhat hi^er stiffness near 
the tip would tend to reduce the severity of the aeroelastlc effects 
greatly, since these effects tend to be quite sensitive to the stiffness 
near the tip. Also the combination of the aerodynamic data used in the 
calculations may not be realized on an actual wing. The lift on the 
wing- tank combination was taken for section B (reference 6 ) with and 
without a tank with sealed g^p, but the lift on the tank was obtained 
for the tank on the model with section A and gap open, because no data 
were available for the lift on the tank on a model with section B and 
gap sealed. Also, the use of fins on the tip tank to overcome its 
inherent pitching moment would tend to reduce the severity of the static 
aeroelastlc phenomena. 

The lateral-control power of the wing- tank combination exhibits two 
interesting features. The aileron reversal speed of the wing with the 
tank is sli^tly hi^er than the antisymmetric divergence speed, and 
the rolling moment due to aileron deflection increases with dynamic pres- 
sure. Figure 7(b) of reference 2 Indicates that the dynamic pressures 
required to reverse the lateral control of an unswept wing is propor- 
tional to the reciprocal of the sum of the moment arms e^^ and 62 - 

If 02 is zero the reversal and divergence speeds coincide, and if 02 
is negative; that is, if the center of pressure due to aileron deflection 
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Is forward of the elastic axis, the reversal speed is hi^er than the 
divergence speed hecause the lift due to aileron deflection tends to 
Increase the single of attack. The aileron reversal speed has no 
physical meaning In such a case. As shown In figure 4 of the present 
paper the assumed value of 62 Is negative over most of the span In 
the case of the wing with the tip tank; In the case of the wing without 
the tip tank. It Is positive at the tip region, which Is Instrumental 
In determining the aeroelastlc characteristics of a w^.ng. For the same 
reason the rolling moment due to aileron deflection increases with 
dynamic pressure in the case of the wing with a tip tank hut decreases 
in the case of the wing without a tip tank. 

No dynamic effects have been considered In the analysis of this 
paper, so that nothing qualitative may be said concerning the flutter 
characteristics of the wing with the tip tank nor its dynamic-response 
characteristics In abrupt maneuver. However, there is reason to believe 
the wing with tip tank may well be subject to unfavorable dynamic 
phenomena for some conditions of fuel mass in the tank at dynamic pres- 
sures even lower than these at which static aeroelastlc phenomena become 
important. 

Certain quasi-steady dynamic phenomena can be estimated by meeuis of 
the semirigid concept outlined in a preceding section, for instance, the 
effect of inertia in a pull-out at constant load factor. As long as the 
center of gravity is ahead of the elastic axis the effect of Inertia is 
to relieve the static aeroelastlc phenomena. Since the inertia forces 
are related to the normal acceleration -vdiich. In turn. Is related to the 
lift, there is a definite relation between the inertia and aerodynamic 
forces. If the assumption is made that the tail and the fuselage carry 
no lift, then the dynamic pressure at dynamic divergence - that is, at 
divergence under conditions "vdilch permit the airplane as a whole to 
accelerate In a direction normal to the fli^t path - can be estimated 
by multiplying the static divergence speed by the factor 

1 

1 _ ^ Mg 

e, "c W 
•^r 

where clj^ Is the distance of the center of gravity of the wing plus tip 
tank ahead of the elastic axis of the wing, M the mass of the wing plus 
that of the two tip tanks, and ¥/g the mass of the airplane including 

that of the wing and of the two tip tanks. ^For ■ — ^ > 0.2 this 

\ ®lr‘= ^ \ 

correction tends to yield values of qj^ which are somewhat too high. J 

Quasi-static dynamic effects can then be included approximately in the 
static aeroelastlc results presented in figures 3 and 7 using the 
value of q*jj corrected in this manner in the ratio q*/q*j) used as 
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a parameter the abscissa, respectively. In these two figures. A 
similar but more complicated correction factor which takes Into account 
the lifts on the tall and fuselage may be devised. 


Swept Wing ¥lth Boomr-Mounted Lifting Surface 

The spanwlse lift distribution corresponding to three angle- of- 
attack conditions of the swept wing with and without two boom-mounted 
lifting- surface configurations are shown in figure 5 subsonic speeds 
and for dimensionless dynamic pressures q.* of 0 and 0.l69. Similarly, 
the spanwlse lift distributions of the swept wing with and without one 
lifting- surface configuration are shown in figure 6 for supersonic 
speeds and for dimensionless dynamic pressures q of 0 and 2 . 17 . In 
both figures 5 6 the lifting surfaces are considered to be mounted 

on a rigid boom and geared to the aileron with a gear ratio of 1. The 
dimensionless dynamic pressures of q,* = O.I 69 and q = 2.17 both 
represent the negatives of the dynamic pressures which would diverge 
the wing without a lifting surface at subsonic and supersonic speeds, 
respectively. The dimensionless dynamic pressure q is used for the 
supersonic case because e^^^ was taken as 0 for that case, so that q* 

is 0 regardless of q. The antisymmetric lift distributions are plotted 
in the form cc^^cC^^, which is similar to the form cc^^cCj^ used for 

the symmetrical cases; the coefficient C-r is the negative of the con- 
n'd 

ventionally defined coefficient of damping in roll. 

As may be expected, the aeroelastlc effect on the spanwlse lift 
distributions is very large at the relatively high dynamic pressures 
represented in figures 5 and 6 . The effect of the boom-mounted lifting 
surfaces, however, is almost negligible except near the wing tip and 
except in the case of the lift distribution due to aileron deflection. 


The lift coefficients, aerodynamic- center locations, rolling- 
moment coefficients, and wing-tip helix angles obtained by integrating 
the lift distributions shown in figures 5 and 6 are represented in fig- 
ures 8 and 9* As indicated in figure 8 for subsonic speeds, the effect 
of the lifting surface on the lift coefficient is negligible up to the 
highest dynamic pressures likely to be of Interest, that is, for values 
of q* between 0.2 and 0.3- The effect of the lifting surface with 

= 1.5 oa "the aerodynamic-center shift is negligible, but the lifting 
^t 

surface with — =2.0 does have a favorable effect on the aerodynamlc- 


center shift; for q* = 0.2 the aerod 3 mamlc- center shift due to aero- 
elastlc action is O.I 7 for the wing without a lifting surface and for 
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the wing with the lifting surface with 
the wing with the lifting surface with 


— = 1.5 but is only 0.l4 for 

= 2 . 0 . 


For the particular sweptback wing under consideration the rolling- 
moment coefficient due to aileron deflection is substantially Increased 
by the boom-mounted lifting surfaces at the highest dynamic pressures : ' 
of Interest. At q* = 0.2, for instance, the rolling-moment coefficient 


is increased about 50 percent by the lifting surface with 


about 100 percent by the lifting surface with 


-^ = 2 . 0 . 
'^t 


^= 1.5 and 
These increases 


are reflected in similar Increases in the wing-tip helix angle per unit 
aileron deflection. 


At dynamic pressures much hi^er than that corresponding to q* = 0.2, 
the wings with lifting surfaces may diverge if the values of q* given in 
table 2 are approached. For the wing without a lifting surface the 
smallest value of q*^^ is negative (q* = - 0 .l 69 j, and the next larger 

one is also negative so that divergence is iii 5 )Osslble. As may be deduced 
from figures 5^ 6, 8, and 9 the divergence of the wings with lifting 
surfaces is a very localized phenomenon,, affecting only the region of the 
wing near the tip. The aileron reversal speed of the wing with lifting 
surfaces tends to be much hl^er than that of the wing without lifting 
surfaces. (See table 2.) 


As shown in figure 9^ the effects of boom-mounted lifting surfaces 
on the aeroelastic behavior of this sweptback wing at supersonic speeds 
are very similar to the effects at subsonic speeds. The effects on the 
lift coefficient and aerodynamic-center shift are very small for the 

lifting surface with ^ but the rolling-moment coefficient and 

the wing-tip helix due to unit aileron deflection are increased con- 
siderably. The divergence speed of the wing with lifting surface is so 
hl^ as to be of no practical interest, but the aileron reversal speed 
is relatively lower, compared to that of the wing without lifting sur- 
face, than in the subsonic case. 


The lateral-control characteristics shown in figures 5 # 6 , 8 , and 9 
are for surfaces geared to the aileron with a 1:1 ratio. When the sur- 
face is not geared to the aileron the lift-curve slope, aerodj-namlc- 
center location, coefficient of damping in roll, and divergence speed 
are the same as when it is geared. The rolling moment and wing-tip helix 
angle due to aileron deflection as well as the reversal speed are even 
lower, however, for the ungeared surface than they are for the wing with- 
ou. a ... at aataoat. apee.a yc. 


6/ ‘'Sr 
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at q* = 0.169 is 0.261, 0.169, 0.338 for the wing -without a 

lifting surface, -with ungeared surface, and -with geared surface respec- 
tively; similarly, at supersonic speeds at q = 2.17 is 


^8/ ‘80 

0.276, 0 . 207 ,_and 0 . 430 , respectively, for these three cases. The vadues 
of q-**- and q for reversal given in table 2 corroborate -this trend. 
Inasmuch as an ungeared surface does not greatly improve the stability 
characteristics (aerodynamic-center shift) and is responsible for a 
deterioration of the lateral- control characteristics it -will not be con- 
sidered any further. In the folio-wing discussion -the lifting surface 
-will be assumed to be geared to the aileron, the gear ratio being 1:1 
not because this is necessarily the optimum value but because that is 
the value for which the calculations described in this paper have been 
made. 


The results presen-ted so far for the swept back wing -with lifting 
surfaces have been for surfaces -with an area ratio = 0. 02 mounted 
on idealized rigid booms. The effects of changes in lifting-surface 
area (or lift-curve slope) and in boom flexibility are shown in figure 10 
for subsonic speeds. This figure shows that in order to decrease the 
aerodynamic- center shift due to aeroelastlc action below that of the -wing 
-wi-fchout a lifting surface a moment-arm ratio dy’c^ of 1. 5 or nrare is 
required regardless of the area of the surface, unless the boom is quite 
flexible. An Increase in the moment-arm ratio from 1.5 to 2 or a decrease 
in the boom stiffness from infinite rigidity to a -value of q*^^ of 

about 0 . 4 serves to decrease the aerodynamic- center shift more -fchan a 
doubling of the surface area (from K^. = 0.02 to = 0 . 04 ). 

Figure 10 also shows -that an incraase of about 40 percent may be 
had in the -wing-tip helix angle due to unit aileron deflection of -fche 

-wing alone by ad.ding a lifting surface -with K_ =0.02 and — = 1.5 

on a rigid boom. By increasing the area ratio to = 0.04 or by 
reducing the rigidity of the boom until is about 0.4 an additional 

60 -percent Increase may be had, but by increasing the moment-arm ratio of 
the lifting surface from ^ = 1.5 io ^ =2.0 only an additional 
30 -percent Increase is ob-balned. 

The dimensionless dynamic pressures required for divergence and 
reversal of the wing--with-llfting-surface combinations represen-ted in 
figure 10 are given in table 2 . For the combinations -with large moment 
arm, surface area, or with very flexible booms, divergence of the local 
type mentioned previously is likely to occur at relatively low dynamic 
pressures, in some Instances so low as to be of practical concern. The 
reversal speed of all configurations is far too hl^ to be of interest. 
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The effectiveness of a hoom-moimted lifting surface as an 
aeroelastic-effect relieving device is probably best illustrated by 
figure 10. For the case considered in figure 10, that is, the svept- 
back wing flying at subsonic speeds with a value of q* of 0 , I 687 , the 
aerodynamic center is shifted 15 percent rearward from the rigid-wing 
position. As shown in figure 10 for a lifting surface with an effective 
area ratio of 0.02, a moment-arm ratio of 1.5^ and a flex- 
ible boom with equal to 1/3 ('ifeich is twice the value of q* 

considered in the figure) this shift is reduced to 10 percent. Larger 
values of and d/c^ and lower values of are likely to be 

impractical because of dynamic (primarily flutter), mechanical, and 
weight considerations. In varying these three lifting-surface param- 
eters it appears that more benefit may be had by varying the moment-arm 
ratio than by varying the area ratio a corresponding amount but that 
unless the moment-arm ratio is larger than about 1.5 no improvement in 
the shift of the aerodynamic center is had at all. A substantial 
improvement in the shift of the aerodynamic center c£ui be obtained by 
Increasing the flexibility of the boom, but too flexible a boom can lead 
to localized divergence of the wing, as well as to divergence of the 
boom proper; as shown in table 2 the wing diverges irfien q* is 0. 2J7 

and 0.217 in the case of the lifting surface with = 0 . 02 , ^ = 1 . 5 » 

and Q.* 2 jj equal to 1/2 and 1/3^ respectively. The use of a flexible 

boom is also likely to Introduce flutter problems . 


Figure 10 indicates that the lateral-control power eind maneuvera- 
bility characteristics may also be in^roved substantially by a geared 
lifting surface; by using a lifting surface with gear ratio K = 1, 


= 0 . 02 , ^= 1 . 5 , and q* 


BD 


= the wing-tip helix angle is twice 


that of the wing without a lifting surface. Again, a variation in the 
moment-arm ratio appears to be more effective than a proportional 
increase in the area ratio but, again, a minimum value of — ^about — 


in this casej is required to obtain any improvement at all. In general, 

the improvement in the lateral-control characteristics obtainable by 
means of a boom-mounted lifting surface appears to be larger than the 
Improvement in the shift of the aerodynamic center. 


In evaluating the results discussed in the preceding paragraphs 
several facts must be kept in mind. Concerning the specific calcu- 
lations described in this paper, as pointed out in connection with the 
calculations for the tip tdnk, the assumed wing stiffnesses may be 
relatively too low near the tip con^jared with actual practice, so that 
the magnitude of the various static aeroelastlc effects may be over- 
estimated somewhat in these calculations. 
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Furthermore, in the calculations for the hooam-mounted lifting sur- 
face the effects of the upwash of the vlng on the lift of the surface 
and of the downwash of the surface on the lift of the wing tip have been 
neglected. In the case of a rigid boom, the effects of the upwash of the 
wing on the lift of the surface can be taken into account by multiplying 
the lift-curve slope of the lifting surface by a factor which is 

one plus the value of the upwash angle per unit angle of attack of tbe 
wing tip. The upwash angle can be calculated by means of the charts of 
reference 8. Similarly, in the case of a rigid boom, the effect of the 
downwash of the lifting surface on the lift of the wing tip can be taken 
into account by calculating a- factor which is equal to one minus 

the downwash caused by the surface on the three- quarter- choitL line of 
the wing at the wing tip; again, the charts of reference 8 can be used 
to calculate this downwash angle, if desired. The elements in the last 
column of the aerodynamic- influence- coefficient matrix are then 

multiplied by this factor. In the case of a flexible boom the method 
of analysis presented in this paper must be modified sll^tly; for 
instance, the angle of attack of the lifting surface is then equaJ. to 
the angle of boom deformation plus the product of the angle of attack 
of the wing tip and the aforementioned factor tj^. 

Finally, no dynamic effects have been taken into account in the 
calculations, nor can a simple correction be given for quasl-statlc 
dynamic effects. However, qualitatively the quasi-static dynamic effects 
are adverse, inasmuch as they decrease the normal force available from 
the lifting siirfate. The essentially dynamic phen^anena, such as encoim- 
tered in flutter, gusts, or abrqpt maneuvers are also likely to be affected 
adversely by boom-mounted lifting surfaces, particularly by heavy sur- 
faces with long or flexible booms. In general, all means of Improving 
static aeroelastlc characteristics by balancing the effects of bending 
and twisting deformations, rather than by stiffening the structure, 
have certain difficulties in common. Exact balance is difficult to 
achieve, and if it is achieved for one Mach number it may not hold at 
others; certainly a condition of balance obtained at subsonic speeds is 
unlikely to carry over to supersonic speeds. Hor does such a means of 
improving static aeroelastlc characteristics necessarily li^rove dynamic 
characteristics; in fact, more often than not, it effects the dynamic 
characteristics adversely. 

As a result of these considerations no optimum boom-mounted lifting- 
surface configuration can be selected. Such a configuration depends on 
the magnitude and nature of the aeroelastlc effects that must be allevi- 
ated and the wei^t penalty that can be tolerated in order to achieve 
this alleviation. Even for a specific case the static calculations 
descidbed in this paper cannot furnish a con 5 )iete answer, because from 
a static point of view a surface with as large an area as possible on a 
boom, as long and flexible as possible without Incurring local divergence 
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would be desirable, whereas from a dynamic point of view these very 
parajneters are those that may have to be avoided. A small area ratio 
is likely to result in a- relatively Ineffective lifting surface, whereas 
a greater ratio is likely to be inefficient, in that the relatively 
small additional alleviation of static aeroelastic effects which it can 
produce is likely to be overshadowed by the severity of the dynamic 
phenomena for which it may be responsible as a result of its greater 
mass and area. Before an qptimiim or compromise configuration can be 
decided upon, several conflgumtions with booms of varying lengths and 
stiffnesses will therefore have to be analyzed for their static and 
dynamic aeroelastic characteristics. 


CONCLUDING REMARKS 


A matrix- integration method has been presented for calculating the 
static aeroelastic characteristics of a wing with concentrated aero- 
dynajnlc forces at its tip due to tanks or bocm-mounted lifting-surfaces. 
A simplified method of calculation applicable to certain cases has also 
been presented, which is based on the concept of the semirigid wing and 
utilizes the characteristics of the wing alone. 

Some static aeroelastic .characteristics have been calculated for an 
unswept wing with a tip tank and for a sweptback wing with several con- 
figurations of boom-mounted lifting surfaces. The results of these 
calculations indicate that a tip tank is likely to affect the static 
aeroelastic characteristics of an unswept wing adversely and that a 
boom-mounted lifting surface geared to the aileron tends to relieve the 
adverse static aeroelastic characteristics of a sweptback wing; the 
shift of the aerodynamic center and particularly the loss of rolling 
speed can be reduced in this manner. In the in^irovement of these char- 
acteristics the length and flexibility of the boom are found to be some- 
what more effective than the area of the lifting surface. The amount of 
relief of adverse static aeroelastic phenomena is likely to be limited 
by dynamic effects Introduced by the use of these Hiding surfaces, but 
no such effects have been taken into account. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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TABLE 2.- DYNAMIC-PEESSUEE PAEAMETEES AT DIVEI?GENCE 
AND AT AHJEBON EEVEESAL FOE SWEPTBACK WING WITH 
BOCM-MOUNTED LIFTING SUEFACES 
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fb) Linear anf/sy/vmafrical fwisf(dampmg-in-rd( easel 
(c) Aileron del lection (oraS=D. 

Figure 3*- Lift distributions due to uniform, angle of attack, linear 
antisymmetric angle of attack, and aileron deflection for unswept 
wing. (M = 0.8.) 
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(a) Unswepf wing. 

(b) 5wepfback wmg. 


Figure 4 .- Spanwise variation of dimensionless moment arms of lifts due 
to angle of attack and due to aileron deflection for swept and \mswept 
wings. 
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(a) Uniiorm angle of aiiack. 

(b) Linear anfisymmefncol fw/st [damping-in-roll case). 

(c) Aileron deflection Ca^6 = l). 


Figure 5-- Lift distributions due to uniform angle of attack, linear 
antisymmetric angle of attack, and aileron deflection for sweptback 
wing (stibsonic speeds, Kv = 0-02). 
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&7j Uniform angle of attack. 

(b) Linear anfisymmefrical fwisf Ldampmg-in-rolJ case). 

(c) Aileron defleciion (a^6=t). 


Figure 6.- Lift distritutioas due to unifonn. angle of attack, linear 
antisjrmmetric angle of attack, and aileron deflection for sweptback 

wing ^supersonic speeds, Ky = 0.02, — = I. 5 V 
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Figure 7*“. Llft-coefflcient ratio, lateral center of pressure, damping- 
ln-roll*coefficlent ratio, rolling-mcment-coefficient ratio, and wing- 
tip helix angle per unit aileron deflection for unswept wing with and 
without tip tank (M = 0.8) . (Note: q* Is based on CLq^ of wing 

without tank.) 
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With boom- mounted lifting surface 

wing alone 


Divergence 



.1 .Z .3 .4 

Dimensionless dynamic pressure, 


Figure 8.- Lift- coefficient ratio, shift in aerodynamic center, rolling- 
moment-coefficient ratio, and wing-tip helix angle for swept wing with 
and vithout geared 'boom-moimted lifting surfaces. (Subsonic speeds, 

0 
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Figure 9 .- Lift-coefficient ratio, shift in aerodynamic center, rolling- 
moment -coefficient ratio, and wing -tip helix angle for swept wing 
with and without geared boom-i^unted lifting surfaces. /Supersonic 


speeds, Kv = 0 . 02 , 


4 ■ 
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Figure 10.- Tlie effect of lifting- surface size and longitudinal location 
on the Ilf t-coefflc lent ratio, aerodynamic-center shift, rolling- 
moment- coefficient ratio, and wing- tip helix angle for a given 
dynamic pressure, (q* = O.I 687 , subsonic speeds.) 


NACA-Langley - 8-27-G2 -3S0 


3EC0R1”1~V' ns4F'C> 



